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(^ ■ ABSTRACT 

o '. 

r^ , Context. Though still beyond our observational abilities. Population III stars are interesting objects in many perspectives. They are 

responsible for the re-ionisation of the inter-galactic medium. They also left their chemical imprint in the early Universe, imprint 
which can be deciphered in the most metal-poor stars in the halo of our Galaxy. 

Aims. Rotation has been shown to play a determinant role at very low metallicity, bringing heavy mass loss where almost none was 
expected. Is this still true when the metallicity strictly equals zero? The aim of our study is to get an answer to this question, and to 
determine how rotation changes the evolution and the chemical signature of the primordial stars. 

Methods. We have calculated seven differentially-rotating stellar models at zero metallicity, with masses between 9 and 200 Mq. For 
(-H ' each mass, we have also calculated a corresponding model without rotation. The evolution has been followed up to the pre-supemova 

P |! stage. 

I ■ Results. We find that Z = models rotate with an internal profile n(r) close to local angular momentum conservation, because of a 

O ' very weak core-envelope coupling. Rotational mixing drives a H-shell boost due to a sudden onset of CNO cycle in the shell. This 

j—j ' boost leads to a high '*N production, which can be as much as 10* times higher than the production of the non-rotating models. 

C/5 , Generally, the rotating models produce much more metals than their non-rotating counterparts. The mass loss is very low, even for 

Cu . the models that reach the critical velocity during the main sequence. It may however have an impact on the chemical enrichment of 

the Universe, because some of the stars are supposed to collapse directly into black holes. They would contribute to the enrichment 

only through their winds. While in that case non-rotating stars would not contribute at all, rotating stars may leave an imprint in their 

►^ , surrounding. Due to the low mass loss and the weak coupling, the core retains a high angular momentum at the end of the evolution. 

J^ . The high rotation rate at death probably leads to a much stronger explosion than previously expected, changing the fate of the models. 

j: ' ' The inclusion of our yields in a chemical evolution model of the Galactic halo predicts log values of N/0, C/O and '^C/'^'C ratios of 

[Zr I -2.2, -0.95 and 50 respectively at log 0/H + 12 = 4.2. 

("^ ' Key words, stars: evolution - stars: rotation ~ stars: chemically peculiar supemovae: general 

o 

OO . 1 . Introduction become strong enough to play an important role in the accretion 

^^ ' , . --- ,_ ---. , .... , rate. This again would reduce substantially the final mass range 

' . Population III (PopIII) stars occupy a key position in the evolu- ^^ ^j^^ ^^.^^ ^^^^ ^j^^^ j^ -^ ^^^ ^^^^ ^^^^ suggested that the most 

-^ ' tion of the Universe. They are the first sources of light that have ^^^^^ .t^,, known to this day may be low-mass PopIII that 

re-ionised the Universe. Also they are the first producers of the j^^^^ undergone surface pollution (Shigeyama et al. 2003), it is 

^ ■ heavy elements in the early Universe'. ^^^ ^-^^^ ^^^^ ^^^ ^^^^ metal-free stars were too massive to 

_Cd Because of the pectihar conditions of the early Universe, es- ^^^^ ^^ ^^^^^^ nowadays (M > O.8M0). The observation facil- 

pecially the lack of coolants provided by the metals, the very first .^.^^ ^^ ^^^ ^^j^ ^^ ^^^^^ ^^ ^^^^^^ ^^^ ^ ^^^^^ ^^^^^^.^^^ 

metal-free stars that forined are supposed to be very inassive, ^^ ^^^^ ^^^ ^^^ ^g^^^^. ^^^ ^^ pair-instability supernovae at 

""'"^ ""l^^TlfT' ^ " i°T r ^° ^®inn? "' "'■ 2002; Broinin intermediate redshifts, see Scannapieco et al. 2005). The only 
et al. 2002). Nakamura & Umemura (2001) propose a bimodal . » • »u ■ u • 1 »• • • j- » u 

,„ ,^ . , , . , ,„->., , , , . , , . , way to constrain their physical properties is indirect observa- 

IMF, with oiie peak at M ^ 10^ Mq and the other at M ^ 1 Mq. ^.^^^ ^^ ^j^^.^ ^j^^^^^^l ^^ ^^^.^^^^ ^^ ^^^^^^ ^^^ ^^ .^ .^ ^^^ 

The bimodality is due to a threshold in the initial density of the ^^^j ^^^^^^^^^^^ whether the most metal-deficient stars of the 

cloud. Since this first siiggestion that PopIII may not all be very j^^j^ ^^ ^^^. ^^^ ^^^ ^^^j ^^^^^^ generation stars (Salvador! 

massive objects, several studies show that the radiation field of ^^ ^j 2007), most of the authors consider they are born from a 

a pre-existing star have a strong impact on subsequent star for- ^j^^^ ^^^.^^^^ ^ ^^^ ^ ^^^ primordial progenitors (Chrisdieb 

mation (Omukai & Yoshii 2003; O Shea et al. 2005; Greif & ^^ ^j 2002; Umeda & Nomoto 2003). Recently, accurate mea- 

Bromm 2006). This would allow lower-mass metal-free stars to ^^^^^^^,^ ^f their surface abundances have started to allow some 

appear Moreover Silk & Langer (2006) show that even if the constraints to be put on the theoretical results of PopIII evolution 

early Universe was free of magnetic fields, magnetic seeds can ^^^ .^j ^^ ^j 2004; Barklem et al. 2005; Cohen et al. 2007; Aoki 

be generated in the disc surrounding an accreting stellar progen- j. , ono?") 

itor. These seeds can be amplified through a dynamo process and '.',.. „ ttt 
o .- r jjj view of their importance, PopIII stars have been the sub- 






Send offprint requests to: Sylvia.Ekstrom@obs.unige.ch ject of many evolutionary calculations. Amongst the latest ones, 

' According to the standard Big Bang Nucleosynthesis, only ele- most are presenting the evolution until AGB phase or central C 

ments with nucleon number A < 7 are formed at a significant level ignition (Chieffi et al. 2001; Marigo et al. 2001, 2003; Siess et al. 

during the first minutes of the Universe (e.g. locco et al. 2007). 2002; Gil-Pons et al. 2007; Lau et al. 2008). Others ai'e focused 
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on final stages (Heger & Woosley 2002; Umeda & Nomoto 
2005). To our knowledge, there are two studies that include ro- 
tating models: the one by Marigo et al. (2003), and the one by 
Heger et al. (2000b). Marigo et al. (2003) study very massive 
objects with masses between 120 and 1000 Mq. In their work, 
rotating stars are treated as rigid bodies and the effects of rota- 
tion intervenes mainly through a correction factor applied on the 
radiative mass loss (see Sect. 2.2). In these models, no rotational 
mixing is accounted for. Heger et al. (2000b) follow the evo- 
lution of differentially-rotating stellar models of 15 - 250 Mq 
from the zero-age main sequence (ZAMS) to core collapse, but 
the detailed results of these models have not been published. 

The present paper explores the effects of differential rota- 
tion on primordial massive stars. We think it is worth doing 
such a study for two main reasons. First, for metallicities be- 
tween 0.004 and 0.040, the inclusion of the effects of differen- 
tial rotation in stellar evolution codes has allowed to improve 
the agreement between the models and observations in many as- 
pects: it has allowed to better reproduce the observed surface 
abundances (Heger & Langer 2000; Meynet & Maeder 2000), 
the ratio of blue-to-red supergiants at low metallicity (Maeder 
& Meynet 2001), the variation with the metallicity of the WR/O 
ratios and of the numbers of Type Ibc to Type II supernovae 
(Meynet & Maeder 2003, 2005). It is interesting to see what a 
code that has been successfully confronted to observations can 
predict in a metallicity domain where no direct observations are 
available. Second, rotation has been shown to have important ef- 
fects on low metallicity stars (Meynet et al. 2006). During the 
Main Sequence (MS), the radiative winds are very low because 
of the lack of metals, so the star cannot get rid of its angular mo- 
mentum. If it starts its evolution with a sufficiently high equa- 
torial velocity, it may reach the critical limit. In that case, the 
outer layers are no more bound and some mass is lost, proba- 
bly through a decretion disc (Owocki 2005) that may be quickly 
erased by the radiation. Further in the evolution, the rotational 
mixing brings freshly synthesised metals at the surface. The 
opacity of the outer layers are enhanced and the radiative mass 
loss is increased. Whereas the low metallicity stars were sup- 
posed to lose almost no mass at all, it has been shown that fast 
rotation can dramatically enhance the mass loss and thus modify 
drastically the evolution of these stars. How does a total absence 
of metals alter the above picture? It is interesting to extend the 
study to strictly Z = models. 

The paper is organized as follow: in Sect. 2, we describe the 
physical ingredients of the modelling; Sect. 3 is dedicated to the 
presentation of the evolution of the models; in Sect. 4, we present 
the pre-SN structure and the yields obtained; we study the effect 
of their inclusion in a galactic chemical evolution code in Sect. 5, 
and then a discussion follows in Sect. 6. 

Let us mention that electronic tables for the evolu- 
tionary sequences and the final structure files are avail- 
able at http://obswww.unige.ch/Recherche/evol/ 
Rotating-primordial- stars. 



2. Physical ingredients of the models 

The stellar models presented here have masses ranging from 9 
to 200 Mq. The initial composition is X = 0.76, Y - 0.24 and 
Z = 0. The models with masses between 25 and 85 Mq were 
computed until the end of the hydrostatic core Si burning. The 
9 Mq model has developed a strongly degenerated core before 
carbon ignition and was thus stopped then. The 15 Mq model 
was stopped at the end of oxygen-burning also because of a too 



Table 1. Rotational properties of the models. The masses are in solar 
mass units, and the velocities in km s"'. The last column gives the cen- 
tral hydrogen mass fraction when the model arrives at the critical limit. 



Mass 


f'mi 


v/v„,i 


u(MS) 


''''c.b 


9 


500 


0.54 


430 


- 


15 


800 


0.72 


626 


0.018 


25 


800 


0.63 


657 


0.253 


40 


800 


0.56 


700 


0.326 


60 


800 


0.52 


740 


0.360 


85 


800 


0.48 


754 


0.399 


200 


800 


0.40 


733 


0.428 



degenerate core at that time. Because of numerical difficulties, 
the 200 Mq was stopped shortly after central He exhaustion. 

The computation was carried out with the Geneva stellar evo- 
lution code. The version of the code is the one used in Hirschi 
et al. (2004) that includes the necessary implementations to go 
on through advanced stages up to the pre-SN. It has been slightly 
modified to join the reaction networks of H and He burning dur- 
ing the MS, as required for zero-metallicity models. 

The reaction rates are taken from NACRE database 
(http : //pntpm . ulb . ac . be/Nacre/barre jdatabase . htm). 
The opacities come from Iglesias & Rogers (1996), com- 
plemented at low temperature by the molecular opacities of 
Alexander & Ferguson (1994). The treatment of convective 
instability is done according to Schwarzschild criterion. As 
in Hirschi (2007), convection is treated as a diffusive process 
from oxygen burning onwards. An overshooting parameter 
ttover - 0.2 Hp is used for H- and He-burning cores only. 

2.1. Rotation 

The effects of rotation are computed in the framework of the 
Roche model (see Ekstrom et al. 2008b) and according to the 
shellular-rotation hypothesis (see Zahn 1992). The code couples 
the centrifugal force to the internal structure equations (Maeder 
& Meynet 2000b). This modifies the hydrostatic equilibrium and 
the surface conditions. The rotation-induced instabilities taken 
into account are the meridional circulation, the secular and dy- 
namical shears. 

- The meridional circulation arises from the impossibility of 
a rotating star to be in hydrostatic and radiative equilibrium 
at the same time (von Zeipel theorem). The vertical com- 
ponent of the meridional circulation velocity is u{r,0) - 
U(r) P2(cos ff) with Piix) the second Legendre polynomial. 
The formulation of the radial amplitude U(r) has been deter- 
mined by Zahn (1992) and Maeder & Zahn (1998) as: 



U(r) = 



pgCpT [V,d - V„d + (<P/SW^] 

I L r , Cpd&\ 



where Cp is the specific heat, Vad - -f^ the adiabatic gra- 
dient, M* - m(i - 2r~)' ^'^ ® ~ P/P i^ ^^^ ^^^^° °f ^^^ 
variation of the density to the average density on an equipo- 
tential. (p and S arise from the generic equation of state: 
— - a^+ip— -6^. En and E/^ are terms which depend on 
the Q- and //-distributions respectively (see Maeder & Zahn 
1998, for details on these expressions). En is the dominant 



Ekstrom et al.: Effects of rotation on the evolution of primordial stars 



term in the big parenthesis, and it can be reasonably approx- 
imated by: iin - | U ~ ^^ i%w)- ^^^ ^^^^ ^" brackets is 
known as the Gratton-Opik term. It may become largely neg- 
ative in the outermost layers of the stars, where p becomes 
very small, which leads to U(r) < 0. In this case, the circula- 
tion goes down along the polar axis from the surface inwards, 
and then rises in the equatorial plane outwards to the surface, 
bringing angular momentum to the surface. 
Differential rotation induces shear turbulence at the inter- 
face of layers having different rotational velocities. A layer 
remains stable if the excess kinetic energy due to the dif- 
ferential rotation does not overcome the buoyancy force. It 
is known as the Richardson criterion and is expressed as 
> Ric, where u is the horizontal velocity, z the 



Ri 



N' 



(dujdzr 



vertical coordinate, and N^ - jj- (Vad - V + (^/5)V^| is the 
Brunt- Vaisala frequency. Ri^ = 1 /4 is the critical value for 
stability. The coefficient of diffusion by dynamical shear im- 
plemented in the code is described in Hirschi et al. (2004): 

Ddyn = -rAQAr 

where r is the mean radius of the zone constituted of consec- 
utive shells with Ri < Ri^ = 1 /4, AQ is the variation of Q 
over this zone and Ar is the extend of the zone. The timescale 
for the dynamical shear is just as it sounds: the dynamical 

timescale Tjyn - -\/rV(Gffv). The prescription for Ddyn is 
valid when the thermal adjustment timescale Tth = -^ (for a 
sphere with diameter d and a thermal difFusivity K) is longer 
than the dynamical timescale. 

Secular shear arises when the thermal dissipation reduces the 
buoyancy force. The coefficient of diffusion by secular shear 
turbulence is determined by Maeder (1997) as: 



D. 



A{K + Dh) 



shear 



[fv,(i + f) + (v.d-v„d)] 

Hp 



?H"^| -(v-v) 



with K 



4acT^ 



^^ 2r the thermal diffusivity, a the fraction of 
the excess energy in the shear that will contribute to mix- 
ing, and (V - V) the difference between the internal non- 
adiabatic gradient and the local gradient. The latter term is 
extremely small and is neglected in the calculation. D/, is the 
coefficient for the horizontal turbulence. It is determined as 
in Zahn (l992):Dh ^ ^ r \2V{r) -aU{r)\, where c,, is a 
constant of order 1, V{r) is the horizontal component of the 
meridional circulation velocity, U{r) its vertical component. 



and in this expression a - 



1 dXriir^O.) 

2 d\nr 



In the zones where Ri < Ri^, the shear mixing coefficient is even- 
tually set to: 

^shear ~ ^dyn "i" ^shear 

The transport of angular momentum inside a star is im- 
plemented following the prescription of Zahn (1992), comple- 
mented by the works of Talon & Zahn (1997) and Maeder & 
Zahn (1998). In the radial direction, it obeys the following equa- 
tion: 

4 ('^"U^^!^! (-'-«) 41 Hf) 



The first term in the right hand side of this equation is the di- 
vergence of the advected flux of angular momentum, while the 
second term is the divergence of the diffused flux. The coeffi- 
cient D is the total diffusion coefficient, taking into account the 
various instabilities which transport angular momentum (con- 
vection, shears). After the end of the MS, the timescale of the 
evolution shortens, so only the diffusive term is kept in the for- 
mulation of the transport of the angular momentum. 

The transport of chemical elements by the meridional cir- 
culation has been shown to be satisfactorily approximated by a 
diffusion process (see Chaboyer & Zahn 1992) throughout the 
whole evolution. 

Since very low- or zero-metallicity stars are more compact 
than their metal-rich counterparts, a similar ratio of the initial 
velocity to the critical velocity^ w/t'crit leads to a higher equato- 
rial velocity. In order to explore the effects of rotation on metal- 
free models, we chose a high initial velocity, 800 km s"'. Note 
that for the 60 Mq model, this initial velocity corresponds to a 
total angular momentum content that is similar to a solar metal- 
licity model of the same mass rotating at 300 km s"' . The 9 Mq 
model was computed with a lower initial velocity (500 km s"') 
because 800 km s"' is overcritical on the ZAMS for this model. 
Table 1 presents the initial equatorial velocities of the models, 
the corresponding u/vcnt ratio and the mean velocity during MS. 

In order to be able to compare our results with other works, 
we computed, for each mass, a corresponding model without ro- 
tation. 



2.2. Mass loss 

Zero-metallicity stars are supposed to evolve at constant mass 
due to a scaling of the wind with metallicity of the form M oc 
{ZlZof-^. However, according to Kudritzki (2002), H and He 
lines are able to remove some mass if the luminosity is suf- 
ficiently high. We thus adopted the Kudritzki (2002) prescrip- 
tion for the models presenting a luminosity log(L/L0) > 6. 
The 85 and 200 Mq models were computed with it from the 
start, while we turned it on during the MS evolution of the 60 
Mq. Since Kudritzki's prescription does not correspond to the 
strictly Z = case, we used the same adaptations as in Marigo 
et al. (2003): 1) the surface metallicity is set Z = max(Zn,in, Zgff), 
where Z^in = 10""* Z0 is the lowest value in the validity range 
of Kudritzki's prescription, and Zetf - \ - X - Y is, the effective 
surface metallicity of the model (allowing for self-enrichment); 
and 2) the effective temperature is set T^s - min(reff , 60'000 K), 
i.e. we did not extrapolate Kudritzki's formula beyond its valid- 
ity domain but constrained T^s ^ 60'000 K for the calculation of 
the mass loss. The models with a mass lower than 60 Mq never 
reach log(L/LQ) > 6, so they were computed without radiative 
mass loss. 

Rotation, by lowering the effective gravity at the surface, 
helps the radiation removing mass from the star. As shown in 
Maeder & Meynet (2000a), the line-driven mass loss is enhanced 
by a factor which depends on the angular velocity such that the 
mass loss in the case of rotation becomes: 



M(0) = Fq, ■ M(0 = 0) where Fq, 



{\-Y)i 



V 






rf 



^ The critical velocity u^^ is reached when the gravitational acceler 
ation is exactly counterbalanced by the centrifugal force. In the frame- 
work of the Roche model, one has: v^, 
polar radius at break-up limit. 



/2 GM 
V 3 '^pb ' 



where ^pb is the 
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Fig. 1. Evolution of Z = models (with rotation: solid lines; without rotation: dotted lines). Left: Hertzsprung Russell diagram. The grey area 
shows the zone of the diagram where He burns in the core of the rotating models. Right: central temperature vs central density diagram. The ignition 
of the different burning stages are given for the rotating models. The grey area indicates the zone of the diagram where the pair-creation instability 
occurs. The dashed grey line shows the locus of the points where the perfect gas pressure is equal to the completely-degenerate non-relativistic 
pressure. 



We see that this expression depends also on the Eddington factor 
r = L/L^dd = kL/(47tcGM) where k is the electron-scattering 
opacity. 

Moreover, there is another situation in which the models can 
lose mass: it is when the external layers of the star reach the crit- 
ical velocity. The over-critical layers are no more bound to the 
star, though it is not yet clear whether they are expelled in the 
interstellar medium or remain in the surrounding, maybe in the 
form of a disc. However, the removal of these layers brings the 
surface back to sub-critical velocity until the evolution of the star 
accelerates it again. In the present calculation, as in Meynet et al. 
(2006), once the model reaches the break-up limit, we adapted 
the mass loss rate to keep the surface just under the critical ve- 
locity. This treatment was applied to all the models except the 9 
Mq, which never reached the critical limit. 



3. Evolution 

3. 1 . Hertzsprung Russell diagram and central conditions 

In the left panel of Fig. 1, we present the Hertzsprung Russell 
diagram (HRD) of our rotating and non-rotating models, and in 
the right panel, their tracks in the log r^ -logpc diagram. We no- 
tice the expected characteristics of rotating models: the ZAMS is 
shifted towards lower effective temperature and luminosity with 
respect to the non-rotating case^. Then, when the evolution pro- 
ceeds, the tracks become more luminous, and the main-sequence 
turn-off is shifted to cooler temperature: the core of the rotating 



' We recall that this shift is due to the sustaining effect of the rota- 
tion: the gravity is counter-balanced both by the gas pressure and the 
centrifugal force in such a way that the star behaves like a lower mass 
one. 



models is refuelled by fresh H brought by the mixing. It thus 
grows, leading to an enhancement of the luminosity. 

The 9 and 15 Mq models are powered only by pp-chains 
when they arrive on the ZAMS and thus they continue their 
initial contraction. It is only when they have produced a suf- 
ficient amount of carbon through 3a reaction (about 10"'^ in 
mass fraction) that the CNO cycle can start, drawing a hook in 
the log Tc - logpc diagram. The onset of the CNO cycle in these 
models can also be seen in the HRD: their tracks evolve towards 
the blue side of the diagram, until the energy provided by the 
CNO cycle stops the contraction and bends the tracks back in 
the usual MS feature. In the rotating 9 Mq, this happens at an 
age of 12.2 Myr (when the central H mass fraction is X^ = 0.439) 
while in the non-rotating one it happens a little earlier, at an age 
of 10.9 Myr (but at a similar burning stage: X^ - 0.439). In the 
case of the non-rotating 15 Mq model, it happens after merely 
1.5 Myr (X^ = 0.695), while it takes 2 Myr (X^ = 0.677) in the 
case of the rotating one. Let us mention that Marigo et al. (2001) 
find that the mass limit for the CNO ignition to occur already on 
the ZAMS is 20 Mq. Our results are consistent with that limit. 

After central H exhaustion, the core He-burning phase 
(CHeB) starts right away: the core was already hot enough to 
burn a little He during the MS so it does not need to contract 
much further. The transition between core H burning and He 
burning is smooth because the model is continuously sustained 
by core nuclear burning. This prevents the models to start a red- 
ward evolution, so they remain in the blue part of the HRD at 
the beginning of CHeB. The H-burning shell is powered only by 
pp-chains and remains radiative. The main part of the luminosity 
is provided by the core (see in Fig. 2 left the case of the rotat- 
ing 15 Mq). Something particular happens to the rotating mod- 
els during the CHeB phase: because of rotational mixing, some 
carbon produced in the core is diffused towards the H-burning 
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Fig. 2. Internal structure of the 15 Mq rotating model at different stages of the core He burning. The solid line is the He-burning energy production 
and the dashed line the H-burning energy production. The dotted line is the fraction of the total luminosity (scale on the right). The grey areas 
show the convective zones. Left: before the CNO shell boost (age: 13.739 Myr; Y^ = 0.4632) - centre: at the moment of the boost (age: 13.744 
Myr; y, = 0.4919) - right: after the boost (age: 13.757 Myr; Y^ = 0.4840). 



shell, allowing a sudden ignition of the CNO cycle in the shell. 
This boost of the shell leads to a retraction of the convective 
core and a decrease of the luminosity. At the same time, it trans- 
forms the quiet radiative H-burning shell into an active convec- 
tive one (see Fig. 2 centre). Later, the model takes the structure 
of a higher metallicity one: convective He-burning core, convec- 
tive H-burning shell and expanding envelope (see Fig. 2 right). 
This boost appears in Fig. 1 {left) as a deep V-shaped feature in 
the rotating tracks of the 9 and 15 Mq, and as a slight depres- 
sion in the 25 and 40 Mq tracks: the luminosity decreases at the 
same time as the T^s cools down. The only non-rotating models 
that present the same feature are the 25 and 85 Mq. In this case, 
the structure is such that the growing convective core reaches 
the bottom of the quiet radiative H-burning shell, triggering the 
same event than in the rotating models. Such a CNO boost leaves 
an imprint in the log T^ - logpc diagram: a second hook can be 
seen in the tracks, which takes the shape of a complex loop in 
the case of the 9 Mq. 

The crossing of the HRD is slow, contrarily to higher metal- 
licity models. As we mentioned above, the nuclear burning is 
continuous inside the core, therefore the Z -Q models do not go 
through a phase of rapid core contraction and its mirror effect of 
fast envelope expansion. On the HRD, we note that the crossing 
of the diagram takes place during the whole CHeB (grey area). 

In the non-rotating set, all the models with mass between 9 
and 60 Mq end their life as blue supergiants (BSG), while the 
85 and 200 Mq become red supergiants (RSG). This compares 
well with the tracks obtained by Marigo et al. (2001). As in their 
study, all the models with masses between 9 and 60 start and fin- 
ish core He burning in the blue part of the HRD. Only the highest 
masses models (M > 70 M0) end their central He combustion 
on their Hayashi line. 

Generally, the rotating tracks end their evolution in a cooler 
part of the HRD than their non-rotating counterparts. The only 
exception is the 85 Mq, which ends its life as a BSG because 
of a blue loop after CHeB. The rotating 15 and 25 Mq end their 
evolution in the blue part of the HRD, without any redward in- 
cursion. The rest of the rotating models become RSG. Let us 
note that none of our models meet the conditions to be Wolf- 
Rayet stars. The 85 Mq ends its life with a log Teff higher than 
4.3, but its surface abundance of hydrogen remains too high. 



Table 2. Final masses, total lifetime, MS and core He-burning dura- 
tions, masses of the He- and CO-core for our models at the pre-SN 
stage. Masses are in Mq, velocities in km s"' and times in Myr. 



Mass 


V\m 


Mfi„ 


Tlifc 


th 


THc 


MHe 


Mco 


9 





9.00 


22.0 


19.50 


1.99 


1.73 


0.46 


9 


500 


9.00 


27.1 


23.90 


2.58 


2.09 


0.63 


15 





15.00 


11.6 


10.60 


0.83 


3.62 


2.92 


15 


800 


14.96 


14.4 


13.20 


0.98 


4.31 


2.71 


25 





25.00 


7.29 


6.64 


0.56 


7.51 


5.39 


25 


800 


24.75 


8.47 


7.88 


0.50 


9.14 


6.30 


40 





40.00 


5.03 


4.60 


0.36 


15.95 


14.52 


40 


800 


37.99 


5.88 


5.41 


0.41 


17.22 


13.53 


60 





59.99 


3.96 


3.60 


0.30 


26.52 


25.11 


60 


800 


57.59 


4.57 


4.20 


0.31 


30.86 


30.58 


85 





84.93 


3.42 


3.05 


0.31 


34.77 


34.50 


85 


800 


74.57 


3.84 


3.52 


0.28 


44.38 


43.92 


200 





199.95 


2.53 


2.27 


0.23 


101.63 


95.69 


200 


800 


183.56 


2.83 


2.56 


0.24 


109.65 


95.97 



3.2. Lifetiives 

Table 2 shows the lifetimes for the MS and CHeB phases. During 
MS, rotation induces an increase of the lifetime from 12.8% for 
the 200 Mq up to 24.5% for the 15 Mq. The models with higher 
initial vjvaii show a more important increase. For the models 
that reach the critical limit, we find an almost linear relation be- 
tween the relative enhancement of the MS lifetime Ath and the 
initial ratio v/vah of the form 



Ath = 



TH(^'/t^crit) - th(0) 

th(0) 



= 0.362 v/va 



(when Th is expressed in Myr) which reproduces the values of 
Ath obtained in our models with an accuracy better than 7%. 

How does our work compare with the other studies? For 
our non-rotating 9 Mq, the MS lifetime is shorter by 12% than 
that found by Marigo et al. (2001), which is consistent with our 
lower overshoot parameter: the convective core of our model is 
about 10% smaller than theirs. Let us mention also that they 
have a higher initial hydrogen abundance (X=0.77). Gil-Pons 
et al. (2007) provide two models of 9 Mq, one computed with- 
out overshoot and the other with the overshoot theory presented 
by Eldridge & Tout (2004). Their calculation of the convective 
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Fig. 3. Models of 20 Mq at various metallicities, with Qini/ficrit = 0.5. Left: internal profile of U(r), where u(r, 9) the vertical component of the 
velocity of the meridional circulation is uir, 0) = U{r) /'2(cos 6) with P2 the second Legendre polynomial. The radius is normalized to the outer 
one. All the models are at the same evolutionary stage, when the central H mass fraction is about 0.40. Right: evolution of the equatorial velocity, 
normalized to the initial velocity. 



boundaries could be approximated by aover = 0.3 //p (P. Gil- 
Pons, priv. comm.). The lifetime of our model is about 21% 
shorter than the model they calculated with this higher over- 
shoot. But it is also 11% shorter than their model calculated 
without any overshoot. In a previous paper (Gil-Pons et al. 
2005), they show the MS evolution of the convective core and 
of the H and He abundances for this model. We notice that the 
transition from pp-c\\am burning to CNO cycle occurs earlier in 
their model (4.4 Myr) than in ours (10 Myr), and in a much more 
abrupt way (see their Fig. 4). The convective core expands sud- 
denly, the H abundance rises in the core and thus allows a longer 
lifetime. The difference in the time at which the CNO cycle ig- 
nites arises most probably from differences in the reaction rates 
and in the equation of state used. 

In the case of our 15 Mq non-rotating model, we find again 
a difference of about 10% with the work of Marigo et al. (2001), 
consistent with the difference found in the 9 Mq models. Our 
lifetime is also 12% shorter than that calculated by Siess et al. 

(2002) for their (non-rotating) 15 Mq, although the overshoot- 
ing parameter they use is the same as ours. Since their paper is 
focused on low- and intermediate-mass stars, we cannot find in 
it any explanation for this difference. 

In order to compare our results with the one of Marigo et al. 

(2003) at higher masses, we computed the evolution of a non- 
rotating 120 Mq up to the end of the MS. Once more the MS 
lifetime of our model is found to be 10% shorter than the one 
computed by the latter authors. 

3.3. Evolution of the equatorial velocity 

Our models are started on the ZAMS with a flat rotation profile. 
As described in Meynet & Maeder (2000), the meridional circu- 
lation establishes then very quickly (1-2% of the MS time) an 
equilibrium profile of Q(r) inside the star with the core spinning 
faster than the surface. 

After this first adjustment phase, the evolution of the sur- 
face (equatorial) velocity results of a delicate interplay between 
the mass loss (which removes angular momentum at the surface) 
and the meridional circulation (which brings angular momentum 
from the core to the surface). How does this interplay behave at 



Z - 0, and what are the differences with the situation at higher 
metallicities? To find an answer, we compared four models of 
20 Mq with Qini/Ocrit = 0.50 at various metallicities, models 
which were computed with the same code as our set of Z = 
models^. In the left panel of Fig. 3, we show the internal pro- 
file of the radial component of the meridional circulation. Let us 
focus on the outer cell, which transports the angular momentum 
outward, when U(r) is negative. The amplitude of the meridional 
circulation is a factor 6 higher in the standard metallicity model 
(Z = 0.020) than in the Z = 0.002 one. This factor amounts to 
25 when we compare with the Z = 10"^ model, and reaches 100 
with the Z = one. This illustrates the effect of the Gratton- 
Opik term (see Sect. 2.1) in the expression of the meridional 
circulation velocity Uir): when the metallicity decreases, stars 
are more compact, so the density increases, and thus U(r) de- 
creases. In the right panel of Fig. 3, we see the resulting evolu- 
tion of the equatorial velocity. At standard metallicity, although 
the amplitude of the meridional circulation is large, the loss of 
angular momentum through the radiative winds has the strongest 
effect, and the equatorial velocity slows down. At low or very 
low Z, the meridional circulation is weak, but the mass loss is 
so diminished that the models are spinning up. In the case of 
Z = strictly, there is no mass loss to remove mass and angu- 
lar momentum, but the meridional circulation is so weak that the 
evolution of 0(r) is very close to local angular momentum con- 
servation, Qr^ = constant: because of the natural inflation of the 
external radius, the surface of the model has to slow down. 

The behaviour described above is valid for all the models 
with Mini < 40 Mq: they are slowing down during the MS evo- 
lution, except at the very end of it when they contract briefly 
before central H exhaustion. Above 40 Mq, the models are spin- 
ning up on the MS. This is not surprising, since the meridional 
circulation is higher in more massive stars (which have a less 
dense envelope), allowing the transport of angular momentum 
that is lacking in lower mass models. 

An important remark has to be made here: during the evo- 
lution, because of the inflation of the radius, the value of the 
critical velocity gets lower. So even a slowing down model may 



These models are taken from Ekstrom et al. (2008b) 
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Table 3. Total mass lost and surface abundances of C, N and O at the end of MS phase and at the end of the evolution for the models that undergo 
mass loss and surface enrichment. Masses are in Mq, velocities in km s"' and abundances in mass fraction. The C abundance is the sum of '^C 
and "C, the N abundance is the abundance in '''N and the O abundance is the sum of "'O, "O and "*0. 



M„ 



Miost 



He 



end of MS 
C 



N 



O 



Mios, 



He 



end of evolution 
C 



N 



O 



9 


500 


- 


0.2426 


4.60e-23 


5.91e-21 


2.35e-22 


- 


0.4279 


5.50e-17 


5.61e-15 


1.85e-16 


15 


800 


0.0306 


0.2503 


1.74e-14 


1.81e-12 


5.90e-14 


0.0418 


0.2646 


6.36e-09 


1.83e-07 


3.38e-09 


25 


800 


0.2252 


0.2517 


4.34e-13 


4.02e-ll 


1.20e-12 


0.2513 


0.2615 


3.53e-10 


9.37e-09 


1.86e-10 


40 


800 


1.0070 


0.2654 


2.24e-12 


1.90e-10 


5.18e-12 


2.0074 


0.3045 


5.17e-07 


1.48e-05 


6.16e-07 


60 





- 


0.2400 


O.OOe+00 


O.OOe+00 


O.OOe+00 


0.0001 


0.2400 


O.OOe+00 


O.OOe+00 


O.OOe+00 


60 


800 


2.3200 


0.2906 


6.84e-12 


5.44e-10 


1.39e-ll 


2.4061 


0.5328 


1.23e-05 


3.54e-04 


3.34e-05 


85 





0.0046 


0.2400 


O.OOe+00 


O.OOe+00 


O.OOe+00 


0.0681 


0.4032 


6.40e-04 


3.88e-03 


1.33e-04 


85 


800 


4.4363 


0.3362 


1.74e-ll 


1.31e-09 


3.19e-ll 


10.4239 


0.6141 


7.29e-06 


1.81e-04 


4.14e-05 


200 





0.0350 


0.2400 


O.OOe+00 


O.OOe+00 


O.OOe+00 


0.0509 


0.2423 


7.93e-13 


6.12e-ll 


1.59e-12 


200 


800 


16.4426 


0.5183 


8.82e-ll 


5.94e-09 


1.23e-10 


20.7609 


0.6685 


7.53e-05 


2.09e-04 


1.02e-04 






A. 

i II 




Fig. 4. Evolution of the fl/Qcnt ratio during the MS. All the models start 
the MS with v^^ = 800 km s"' , except the 9 Mq which starts the MS at 
500 km s-'. 



reach the critical velocity at a time of its evolution. This is in- 
deed what happens with our models. In Fig. 4 we present the 
evolution of the n/Ocrit ratio during the MS. As mentioned in 
Sect. 2.1, the 9 Mq model is the only one that does not reach 
the critical limit. For the other models, the higher the mass the 
sooner the critical limit is attained. The last column of Table 1 
gives the central mass fraction of hydrogen when each model 
reaches the critical velocity. 

Once at critical limit, all the models remain at the limit. In 
Fig. 4, the 85 and especially the 200 Mq seem to depart from 
Q/Qcrit = 1, but this is because the limit shown here is only the 
Q-limit, where the centrifugal force alone is taken into account 
to counterbalance the gravity. In the two above models however, 
the radiative acceleration is strong and they reach the so-called 
QF-limit, that is the second root of the equation giving the crit- 
ical velocity: geff [1 - T] =0 (see Maeder & Meynet 2000a, for 
the detailed derivation). The true critical velocity is lowered by 
the radiative acceleration, and though the plotted Q./Q.ait ratio 
becomes lower than 1, these models are actually at the critical 



limit and remain at this limit till the end of the MS. At central 
H exhaustion, the moderate inflation of the radius brings the sur- 
face back to sub-critical velocities. 



3.4. Mass loss and surface abundances 

In Meynet et al. (2006), we show how rotating low-metallicity 
stars may lose more mass than expected mainly through two pro- 
cesses: first the reaching of the break-up velocity, which leads to 
a mechanical mass loss of moderate amplitude through a decre- 
tion disc; second the enrichment of the surface because of rota- 
tional mixing, which leads to a very large enhancement of the 
radiative mass loss. Now what happens if the metallicity is as 
low as strictly zero? The total mass lost at the end of MS re- 
mains modest in the present models (Table 3): between 0.03 Mq 
for the 15 Mq model (which reaches the critical limit only at the 
very end of the MS, when X^ = 0.02), and 16.4 Mq for the 200 
Mq model (which reaches the break-up limit at Xc = 0.43). The 
mechanical mass loss removes only the outermost layers, which 
have a very low density, so only a little mass is removed during 
this phase. 

After the MS, as mentioned in Sect. 3.1, the models remain 
in the blue part of the HRD most of the CHeB phase and then 
their redward incursion is moderate. The outer convective zone 
that develops remains very thin, so the dredge-up does not dig 
deeply enough (see the Kippenhahn diagrams in Fig. 5). Most 
of the post-MS evolution occurs with a surface metallicity much 
lower than 10"^. As a result, the mass lost during CHeB remains 
very low: it amounts to at most 6 Mq in the case of the 85 Mq, 
which is the model that loses the most mass at this stage. 

The total mass lost by the models during the whole evolu- 
tion is small. The 200 Mq model loses 20.7 Mq, the 85 Mq 
model 10.4 Mq. The models have to be as massive as 40 Mq 
or above to lose more than 1 Mq. This picture is very different 
from the evolution with just a small fraction of metals: the 60 
Mq dHZ- 10"** presented in Meynet et al. (2006) loses 36 Mq, 
among which 27 Mq are due to the enhancement of the radia- 
tive winds by the surface self-enrichment. At Z = 0, a similar 
model loses 2.4 Mq only. To check if the choice of a higher ini- 
tial equatorial velocity could enhance significantly the mass loss 
in Z = models, we have computed a test model of 40 Mq with 
yi„i - 900 km s"' {i.e. Sl/Qcrit = 0.83). We have evolved it until 
the end of CHeB. During the MS evolution, it reaches the criti- 
cal limit when its central H mass fraction is X^ - 0.43 and loses 
1.702 Mq until the end of the MS. The total mass lost at the 
end of CHeB amounts to 1 .734 Mq which compares well with 
the 1.571 Mq lost by the model with a slower uini at the same 
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Stage. We can thus conclude that the choice of a higher Uini is not 
determinant to really change the situation about mass loss. 

The non-rotating models do not lose mass at all: even the 200 
Mq model loses only a few hundredths of solar mass. We will 
see in Sect. 4.3 that some of our models may be able to contribute 
to the chemical enrichment of the medium only through their 
wind ejecta. In that context, while the non-rotating models do 
not contribute at all, the rotating ones will leave a modest imprint 
on their surrounding. 

We cannot compare directly our results with the literature: 
the only study that does not compute the evolution at constant 
mass is the one by Marigo et al. (2003) and their masses do 
not match exactly ours. Nevertheless, we can make a qualitative 
comparison. Our rotating 200 Mq model experiences a much 
stronger mass loss than what they found for their rotating 250 
Mq model, which loses 3.35 Mq only, that is 6 times less than 
our slightly less heavy model. Their models have been computed 
with a lower initial velocity, and despite the maximal coupling 
brought by the rigid body rotation, they attain the critical limit 
only at the very end of the MS phase. The mass loss they adopt 
then (Merit = lO"-' M0yr"') is about 2 orders of magnitude 
higher than that resulting from our treatment of the critical mass 
loss, but the period during which they apply it is short. 



4. Final structure and light elements yields 

4. 1 . Kippenhahn diagrams 

In Fig. 5 we present the Kippenhahn diagrams (evolution of the 
structure) of the models between 15 and 85 Mq. Generally, rota- 
tion suppresses the convective instability outside the H-burning 
core by smoothing the abundance gradient of helium above the 
core. The CNO shell boost is clearly visible in the diagrams of 
the non-rotating 25 Mq and the rotating 15, 25 and 40 Mq mod- 
els: the H-burning shell becomes suddenly convective and the 
core is abruptly reduced. We can discover in the rotating 60 Mq 
diagram a similar boost, but only at the beginning of shell He 
burning. This boost left no signature in the HRD. 

As mentioned above, the rotating 15 and 25 Mq models re- 
main in the blue part of the HRD and do not develop any outer 
convective zone (OCZ). The rotating 40 Mq presents a very thin 
one and the 85 Mq a very early and short-lived one. On the con- 
trary, the 60 Mq develops a deeper OCZ at the end of CHeB 
and for almost the rest of its evolution until the pre-SN stage. Its 
surface enrichment is thus larger than that presented by the 85 
Mq model. The only non-rotating model that develops an OCZ 
is the 85 Mq. Though this convective zone disappears quickly, it 
suffices to dredge up heavy elements brought near the surface by 
the CNO shell boost, and it is the model that shows the highest 
level of enrichment at the end of its evolution (see Table 3). 

Limongi et al. (2000) show that at standard metallicity, the 
transition from convective to radiative core C burning occurs be- 
tween 15 and 20 Mq. From the Fig. 1 of Limongi et al. (2001), 
we deduce that at Z = 0, this transition occurs at slightly higher 
masses, between 20 and 25 Mq. We note that our 15 Mq mod- 
els have convective C-burning cores, while cores are radiative in 
the case of our 25 Mq models, in agreement with Limongi et al. 
(2001). In the non-rotating models with M > 40 Mq, the shell C 
burning is convective up to the end of the hydrostatic evolution, 
while in the rotating models, the shell experiences a few flashes 
but otherwise remains radiative. 



H 10- 



- 


■^ 


.^-_-r'-'^ " 


1 / 
1 / ^- ■' ■ '^ 

r no rotation 


- 
rotation 


- 


7, 


'He ' 

12.13C . 


14(4 


'"■'^O 

^^Xe - 



0.1 
0.01 



10 ZO 30 40 50 80 10 20 30 40 50 80 

Fig. 6. Total stellar yields as a function of initial mass for our Z = 
models. Left: models with rotation; right: non-rotating models. 



4.2. Stellar yields 

The yields calculated from stellar models are strongly correlated 
to the value of the CO-core mass, Mqq. We give the values of 
Mco of our models in Table 2. They were determined as in Heger 
et al. (2000a) as the mass coordinate at which the abundance of 
He drops below 10"^ (in mass fraction). From these values, we 
determined the remnant masses of our models, Mrem, with the 
relation between Mrem and Mco described in Maeder (1992) 

In Table 4 we present the total stellar yields for the light 
elements that are not significantly modified by the explosion 
(Limongi & Chieffi 2002). The total yields are the sum of the 
pre-SN yields (mp^''^ ) and the winds yields (m/? "'"'''*). The pre- 
SN contribution from a star of initial mass m to the stellar yield 
of an element / is: 



mp 



preSN 



[Xi(m,) - Xf]dm, 

•J M„m 



with X'^ the initial abundance (in mass fraction) of element / 
(X!' - for all elements except He), and Xi(mr) its final abun- 
dance at the Lagrangian coordinate m^. The winds contribution 

is: 



mPi, 



„ winds 



■T(m) 



M(m,t)[xfim,t)-X^^]dt 



where t(ot) is the final age of the star with initial mass m, M(m, t) 
the mass loss rate when the age of the star is f, and X^ (m, t) the 
surface abundance (in mass fraction) of element / at age t. Most 
of the time, the total stellar yields are dominated by the pre-SN 
contribution. 

We did not include the 200 Mq models in Table 4 because 
they have not undergone the shell He-burning phase, which mod- 
ifies noticeably the results. The 9 Mq models were also stopped 
before central C burning, but as we mentioned, the core was 
strongly degenerated at carbon ignition, which could lead to the 
complete destruction of the star without further processing. Note 
that the carbon ignition occurred at the edge of the core, so there 
is a possibility that the star survived C burning. In that case, since 
the mass included below the He-burning shell is smaller than the 
limiting mass for Ne burning (1.37 Mq), the star is expected to 
collapse in a Type II SN by electrons capture on a degenerate 
O-Ne-Mg core (Nomoto 1984). This would modify the pre-SN 
yields presented here since the C burning would have processed 
the matter of the core. We would thus expect a reduction of '^C 
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Table 4. Total stellar yields at the pre-supernova stage, in Mq. In parenthesis, the "wind-only" contribution is given for the models for which no 
SN contribution is expected (see Sect. 4.3). The masses of the models are in Mq and the velocities in km s"' 



Mini 


Vm, 


in-4 


^He 


'^C 


"c 


14N 


160 


"0 


ISQ 


22Ne 


Z 


9 





-0.53 


1.17 


0.19 


5.54e-ll 


2.68e-8 


0.24 


2.74e-ll 


4.45e-9 


1.58e-7 


0.43 


9 


500 


-0.65 


1.59 


0.24 


3.27e-5 


1.81e-3 


0.17 


2.20e-7 


8.20e-6 


8.98e-6 


0.41 


15 





-1.02 


1.24 


0.22 


5.94e-8 


5.13e-6 


0.86 


1.93e-10 


2.68e-9 


8.73e-7 


2.53 


15 


800 


-1.23 


2.10 


0.61 


9.21e-4 


4.99e-3 


1.05 


8.88e-7 


7.17e-6 


5.44e-4 


2.10 


25 





-1.97 


2.32 


1.26 


2.09e-4 


5.69e-3 


2.78 


1.86e-6 


1.79e-9 


3.27e-6 


4.93 


25 


800 


-2.13 


2.57 


1.32 


3.21e-5 


1.08e-3 


3.76 


2.75e-7 


7.78e-6 


5.54e-4 


5.75 






(-0.02) 


(0.003) 


(5e-14) 


(1.59e-14) 


(6.86e-12) 


(2e-13) 


(1.51e-15) 


(6.51e-19) 


(l.Ole-26) 


(7e-12) 


40 





-3.21 


2.94 


0.23 


4.85e-8 


3.09e-7 


7.25 


2.29e-9 


1.55e-8 


3.63e-6 


10.2 






(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


40 


800 


-3.32 


5.32 


1.63 


8.41e-4 


l.OOe-1 


7.57 


4.04e-6 


6.74e-5 


3.30e-2 


11.4 






(-0.17) 


(0.05) 


(3e-8) 


(7.96e-9) 


(1.07e-6) 


(4e-08) 


(5.89e-ll) 


(5.43e-14) 


(2.64e-ll) 


(le-6) 


60 





-4.87 


4.58 


0.18 


9.07e-7 


2.58e-6 


13.6 


1.69e-8 


2.26e-7 


1.89e-6 


17.7 






(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


60 


800 


-4.76 


3.33 


0.79 


2.39e-4 


1.20e-2 


16.9 


1.89e-5 


1.33e-4 


5.24e-5 


21.7 






(-0.21) 


(0.04) 


(2e-8) 


(7.10e-9) 


(9.00e-7) 


(9e-8) 


(1.94e-10) 


(8.18e-14) 


(4.12e-ll) 


(le-6) 


85 





-7.22 


11.70 


0.28 


1.46e-2 


9.87e-l 


19.7 


6.63e-5 


1.17e-7 


1.18e-3 


25.5 






(-0.004) 


(0.01) 


(3e-5) 


(8.20e-6) 


(2.42e-4) 


(8e-6) 


(7.26e-9) 


(1.50e-ll) 


(3.99e-7) 


(3e-4) 


85 


800 


-6.15 


5.31 


1.00 


6.44e-3 


1.40e-l 


26.2 


1.24e-4 


2.32e-7 


6.47e-6 


31.6 



and an enhancement of '^C, '"*N, '*0 and ^^Ne. In the case of 
the 15 Mq models, the evolution has been stopped at the end 
of central O burning. At that time, the mass interior to the C- 
burning shell is higher than Chandrasekhar mass, so the star is 
expected to undergo an off-centre Si flash and to end with an Fe 
core (Nomoto 1984). We then expect the yields presented here 
to be modified slightly (reduction of C, N and Ne, enhancement 
of O). We have checked the amplitude of the variation when the 
yields are calculated at the end of O burning or at the end of Si 
burning in the models that have been followed until the pre-SN 
phase, so we do not expect the yields of the 15 Mq models given 
in Table 4 to vary by more than 15%. 

Rotation leads to larger yields for almost all the isotopes and 
for almost all the masses (see Fig. 6). The exceptions are the 
non-rotating models that undergo a strong CNO shell boost: the 
25 and 85 Mq, which produce larger amounts of ^^C and ^"*N. 
The 25 Mq overtakes its rotating counterpart also in the produc- 
tion of '^O, and the 85 Mq produces also more "*He and ^^Ne 
than the rotating one. Otherwise, the enhancement of the yields 
brought by the rotational mixing may reach orders of magnitude, 
especially so for the production of primary ^"*N: a factor 10^ be- 
tween the rotating and non-rotating 9 Mq models and even a 
factor 10'' in the case of the 40 Mq models. We note that rota- 
tion makes the production of '"^N robust through the whole mass 
range, while there are enormous differences from mass to mass 
in the non-rotating models. 

Figure 7 shows a comparison of the CNO yields of all our 
Z = models with the yields of the non-rotating models at Z = 
from Chiefs & Limongi (2004) and those of the rotating models 
at Z = 10"** from Hirschi (2007). For '^c (left) and ""O (right), 
our non-rotating 15 and 25 Mq yields compare within a factor 
2 or 3 with those of Chiefli & Limongi (2004). Their 25 Mq 
produces also a huge amount of '"*N (centre), even a factor ~ 10 
higher than ours, so the CNO boost of the H-burning shell seems 
to be a robust feature in this mass domain. Compared to the '^C 
produced at Z = 10 **, the production of our rotating 60 and 85 
Mq models is smaller, while it is almost equal or larger for the 9 
to 40 Mq models. This can be understood because our high mass 
models do not undergo a shell boost, so their He-burning core is 
larger, which leads to a lower '^C yield. The ""O production is 
about the same at Z = 10 ** and Z = 0. Our PopIII models have 



slightly larger '^O yields than the models of Hirschi (2007), but 
we use a higher overshooting parameter which leads to larger 
CO cores and thus larger '*'0 production. 

An interesting point is the production of primary nitrogen 
at very low metallicity. Meynet & Maeder (2002) show that this 
production increases when the metallicity decreases, and Hirschi 
(2007) confirms the trend. We should thus expect a still higher 
production at Z = 0, but this is not the case. '"'N is systemati- 
cally less produced at Z = than at Z - 10"^. We relate this 
change of behaviour to the fact that the Z = stars burn their 
hydrogen at temperatures hot enough to burn some He as well 
and do not evolve into the red part of the HR diagram (or do so 
very slowly during CHeB) as can be seen in Fig. 1 (left). Due to 
this fact, the stars undergo much less mixing during the CHeB 
phase. Indeed, no redwards evolution means that the gradient of 
the angular velocity at the border of the He-burning core remains 
much shallower and triggers much less mixing. On the contrary, 
when a small amount of metals is present, stars begin to burn 
their hydrogen through the CNO cycle from the start. Thus the 
central temperatures during the MS phases are well below those 
required for He-burning. At the end of the MS phase, the core 
contracts and the stars evolve quickly to the red. A steep angular 
velocity gradient appears at the core border and triggers strong 
mixing. Figure 7 shows that there is even a difference of almost 
two orders of magnitude between the production of the 20 Mq at 
Z = 10"^ and the 25 Mq atZ - Q. This can be explained as fol- 
lows: the Z = 10 ** model experiences a second CNO shell boost 
during shell He burning, which enhances its '"^N production by 
about two orders of magnitude. 

4.3. Fate of the models 

At high metallicity, massive stars experience a very strong mass 
loss. Their final mass is therefore very small (even smaller than 
the final mass of stars with an initial mass Mini < 25 Mq, see 
Maeder 1992). Thus, the fate of all massive stars is a SN explo- 
sion, whatever their initial mass. 

At very low metallicity, this is no longer true. Heger & 
Woosley (2002) and Heger et al. (2003) determine the limiting 
Mue for which the fate will be: 

- a Type II SN: Mhe <9Mq, 
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- a black hole (BH) by fallback: 9 Mq < Mee < 15 Mq, 

- a direct BH: \5 Mq < Mhe < 40 Mq, 

- a pulsational pair-instability followed by a SN with BH for- 
mation: 40 Mq < Mhe < 64 Mq, 

- a pair-instability SN (PISN): 64 Mq < Mhc < 133 Mq, 

- or a direct BH collapse: Mhe > 133 Mq. 

Note that the limits above are obtained by non-rotating models, 
and they may not be appropriate to rotating models, as shown 
by Fryer et al. (2001). But since we are lacking more accurate 
determinations, we applied them to all our models. A star that 
ends its life as a BH is supposed to contribute to the chemical 
enrichment of its surrounding only through its winds. We give 
in parenthesis in Table 4 this "wind-only" contribution for the 
models that, according to the limits determined above, end their 
life as BH by fallback or by direct collapse. 

The winds yields are very small, since we have seen that 
the winds are very weak (Sect. 3.4). The 25 Mq for instance, 
releases a total metals amount of less than 10"^ Mq through 
its winds while it produces a little more than 5 Mq in the SN 
scenario. The 60 Mq contributes to the enrichment of the early 
Universe with a mere 10"^ Mq of metals. Note that the net en- 
richment depends of course on the amount of metals ejected, but 
also on the volume in which they are diluted. The wind con- 
tribution is expected to be a slow one, so it will be diluted in a 
smaller volume. The net enrichment may not be so much smaller 
compared to a SN contribution diluted in a larger volume. 

Rotation changes the mass of the He-cores at the end of 
the evolution. They are usually more massive because of the 
rotationally-induced mixing, and this may have consequences 
for the way the star ends its life. While the 25 Mq without rota- 
tion ends as a Type II SN, the rotating one is supposed to produce 
a weak explosion and a BH by fallback. The non-rotating 85 Mq 
remains in the direct BH domain while the rotating one reaches 
a He-core mass that implies pulsational pair-instability, followed 
by a Fe-core collapse into BH. Its true yields must thus be very 
close to the pre-SN value. 



Table 5. Some rotational properties of the iron core of our models: the 
mass Mpe where the abundance of the sum "'**Cr + ''^Fe + '''Ni amounts 
to 50% in mass fraction; the angular velocity Q. taken in the middle of 
the iron core (the profile of f2(m) is almost flat in the core, so this value 
represents the mean value of Q. in the iron core); the maximum specific 
angular momentum j inside the core, corresponding to the first peak 
seen in Fig. 8 {left); the total angular momentum £. and the parameter 
/? = E,otl\Epoi\, both taken at the border of the iron core. 



Mass 


Mpe 


a 


J 


£ 


yS 


[Mq 


[Mq 


[s-'] 


[cm^ s-'] 


[g cm^ s"'] 


% 


25 


1.327 


0.4232 


4.073e+16 


5.837e+49 


1.5833 


40 


1.555 


0.3547 


3.667e+16 


5.963e+49 


1.0046 


60 


2.516 


0.1578 


4.338e+16 


1.152e+50 


0.5449 


85 


2.779 


0.1062 


3.913e+16 


1.182e+50 


0.3624 



As mentioned in Sect. 3.1, none of our models meet the con- 
ditions for becoming WR stars. They end their evolution hav- 
ing kept their envelope, and thus seem to fail to become a GRB 
progenitor as defined in the collapsar model of Woosley (1993). 
However, we have seen that the rotating 85 Mq is likely to un- 
dergo pulsational pair instability^ and thus lose some mass in this 
process. If the mass lost is large, it may become a GRB progen- 
itor since the very low core-envelope coupling has maintained a 
high angular momentum in the core. Fig. 8 (left) shows that in 
the innermost parts of this model, the specific angular momen- 
tum is higher than that needed to create an accretion disc around 
aBH. 

More generally, we expect that rotation will change the way 
a star collapses. Shimizu et al. (2001) show that an anisotropy of 
only 10% in the neutrinos flux can boost the energy of the explo- 
sion. Walder et al. (2005) and Yamasaki & Yamada (2005) show 
that the rotation-induced asphericity of the matter accretion rate 
is an even more important factor that helps producing a robust 
explosion. The geometry of the system leads most probably to 
jet-like explosion. Table 5 gives the rotational characteristics of 



This process has not been followed in the present models. 
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momentum necessary for a non-rotating and a maximally-rotating black hole, respectively. Right: profile of fin, inside the same model. 



the iron core of our models and Fig. 8 (fight) the distribution 
of angular velocity inside the 85 Mq model. All models have 
kept a high total angular momentum, and meet the conditions to 
produce a maximally rotating BH. Though detailed calculations 
would be needed to ascertain this point, the true fate of our mod- 
els could be much more explosive than previously supposed. In 
the chemical evolution model discussed below (see Sect. 5), we 
shall use the total yields, implicitly assuming that the material 
above Mrem will be ejected. 

The centrifugal force has a stabilising effect on the collaps- 
ing core, so the core bounce may occur at sub-nuclear densities. 
This alters the emission of gravitational waves (GW). Ott et al. 
(2004) show that rotation increases the amplitude of GW be- 
cause of the induced anisotropy. But if rotation is really fast, for 
j8 - iirot/l^^potl > 0.3%, the amplitude of GW becomes smaller, 
and the characteristic frequency is shifted to lower frequencies. 
Except for the 85 Mq, our models have a high /3 so their GW 
emission is expected to be damped by their high rotation rate. 
Suwa et al. (2007) show that very massive collapsing PopIII 
stars may contribute to the GW background and be detectable 
by future facilities as DECIGO, the successor of LIGO. In that 
framework, our 85 Mq, with/? = 0.3%, is expected to contribute 
to this background. 

5. Chemical evolution of the Galaxy with Z = 
yields 

Yields of extremely low-metallicity models are constrained by 
the most metal-poor stars observed in the halo of our Galaxy. 
They are supposed to trace the local chemical evolution of their 
surrounding medium at the time of their birth. The notion of time 
is important to keep in mind because stars with different masses 
have different lifetimes: the elements they produce are not re- 
leased at the same moment. The oldest stars observed nowadays 



may have form when only massive stars had time to contribute 
to the enrichment of the medium. If we want to understand the 
abundances observed at the surface of these old stars, we can- 
not just integrate the yields obtained by a set of models with an 
IMF, because we would lose track of the time. What we need is 
a complete chemical evolution model (CEM). 

Such a CEM is computed by Chiappini et al. (2006b) for 
example^. They compare their predictions for the evolution of 
N/O and C/O upon the adoption of different sets of stellar yields 
with the surface ratios determined for the metal-poor stars ob- 
served by Israelian et al. (2004) and Spite et al. (2005). They 
show that the N/O ratio is much higher than what was predicted 
by a CEM using the yields of the slow -rotating Z = 10"^ models 
from Meynet & Maeder (2002) down to Z = 0. When adding the 
yields of the fast-rotating Z = 10"** models from Hirschi (2007), 
the fit is much improved. The same improvement is found for 
the C/O ratio, which presents an upturn at low metallicity. Note 
that the high N/O and the C/O upturn of the low-metallicity stars 
is also observed in low-metallicity DLAs (Pettini et al. 2008). 

Their CEM is computed with the assumption that yields from 
Z = 10"^ stellar models are valid down to Z = 0. But actually, 
there is a physical transition at Ztrans - 10"'": the stars below 
Ztian,s present a PopIII-like behaviour (with some He burning al- 
ready during the MS), and the stars above Ztrans have well sepa- 
rated H- and He-burning phases. 

We computed a similar CEM, but this time including our 
Z = pre-SN yields up to Ztrans, the yields of Hirschi (2007) 
up to Z = 10"^, and then the yields from Meynet & Maeder 
(2002) and Hirschi et al. (2004). For the PopIII stellar models, 
we used an IMF that is truncated at 1 Mq but otherwise stan- 



* The details of the models can be found in Chiappini et al. (2006a), 
where they show that such a model reproduces nicely the metallicity 
distribution of the Galactic halo stars. This means that the timescale for 
the chemical enrichment in the halo is well constrained. 
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dard. Compared to the chemical evolution computed using the 
Z - 10"^ yields down to Z = (Fig. 9, left panel), the new 
computation predicts a reduction of the N/O and C/O ratios at 
extremely low metallicity (Z < 6 ■ 10"^). This trend is expected 
given the higher O and the lower C and N production of the 
PopIII models (Sect. 4.2 and Fig. 7). Note that a similar CEM 
computed with the "wind-only" contribution for our 25 to 60 
Mq models leads to such a slow O enrichment that it would de- 
viate from the CEM with Z = 10"^ yields curve only at a O/H 
domain much lower than that shown in Fig. 9. 

Another interesting quantity is the ratio '^C/'^C. A large N 
production comes with a high '^C value since they are both 
the products of '^C conversion by the CN-cycle in H-burning 
zones. We thus expect a low '^C/'^C ratio from the yields of 
fast-rotating models. Chiappini et al. (2008) use the same set of 
yields as in Chiappini et al. (2006b), but this time to study the 
'^C/'^C ratio. They show that this ratio is indeed much closer to 
the recent observations provided by Spite et al. (2006) upon the 
inclusion of fast-rotating Z = 10"** models. With the inclusion of 
our pre-SN PopIII yields (Fig. 9, right panel), we predict that the 
^^C/'-'C ratio should rise again slightly with decreasing [Fe/H] at 
extremely low metallicity. 



6. Discussion and conclusion 

We computed a set of metal-free stellar models with masses be- 
tween 9 and 200 Mq, following their evolution until the pre- 
supernova stage. The set was composed of seven differentially 
rotating models with initial velocity 800 km s"' (except for the 
9 Mq model which started the main sequence with 500 km s"'), 
and a comparison set of non-rotating models. We showed that in 
PopIII stars, although the effects of rotation are important, they 



are less strong than in stars with even a tiny amount of metals 
(Z = 10"^). There are two main causes for this: 

1 . the meridional circulation is extremely weak at Z = 0. The 
dependance of the outer cell in 1/p (see Sect. 3.3) makes it a 
factor 100 lower compared to the value at standard metallic- 
ity, and still a factor 4 compared to Z = 10"^. The evolution 
of the angular momentum is close to local conservation, so 
the critical limit is reached late. Moreover, only a very small 
amount of mass needs to be lost in order to bring the surface 
back to sub-critical rotation; 

2. there is little structural readjustment at the end of the MS, be- 
cause the core is akeady hot enough to burn some He during 
the main sequence, so the nuclear burnings are continuous. 
The transition between H and He burning is smooth and the 
models remain on the blue side of the Hertzsprung-Russell 
diagram for a long part of the core He burning. When they 
eventually reach the red side, their outer convective zone re- 
mains thin, so the dredge-up is not efficient and the surface 
enrichment stays low. The envelope remains transparent to 
radiation and the mass loss is negligible. 

There is actually a mass loss where none was expected (because 
of the wind dependance on Z), but it remains modest even for 
the most massive of our models. 

Could this picture be modified by changing the physical in- 
puts in the models? We have shown that the situation is qual- 
itatively the same when we take a faster initial rotation rate. 
The main problem is the weakness of the meridional circula- 
tion, which does not efficiently feed the outer layers with an- 
gular momentum taken from the fast spinning core. So, even if 
the critical limit is reached a little earlier, the models still do not 
lose a lot of mass. For this work, we did not take into account 
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the anisotropy of the winds in our models. We do not think it 
would change the picture presented here, because the loss of an- 
gular momentum by radiative mass loss is negligible. There is 
still an interesting possibility: the implementation of magnetic 
fields. They have been shown to change many features in stel- 
lar evolution (Heger et al. 2005; Maeder & Meynet 2005), and 
to provide the strong core-envelope coupling that is lacking in 
the Z = models. With such a coupling, the mass loss becomes 
much larger and the evolution may be strongly altered (Ekstrom 
et al. 2008a). 

With the initial velocities considered here, our models do 
not follow a chemically homogeneous evolution (Maeder 1987; 
Heger & Langer 2000). In the frame of the physics we used, ho- 
mogeneous evolution for Pop 111 stars would need faster rotation 
rates to occur. When the magnetic dynamo mechanism proposed 
by Spruit (1999, 2002) is taken into account, chemical mixing 
becomes more efficient and homogeneous evolution can be ob- 
tained for smaller initial rotation rates (Yoon & Langer 2005; 
Yoon et al. 2006; Woosley & Heger 2006). Note however that for 
extremely low metallicities, the weakening of the meridional cir- 
culation becomes a real handicap for the models to evolve chem- 
ically homogeneously. Hirschi et al. (2008) find that at Z = 10"*^, 
a 40 Mq model (computed with magnetic fields) starting its evo- 
lution with f/vcrit = 0.55 does not rotate fast enough to follow 
such an homogeneous evolution. We expect that at Z = the 
problem would be still more acute. The consequences of such 
extreme effects of rotation is left for a future work. 

It is important to underline here that the only prescription 
for stellar dynamo existing to date, the Tayler-Spruit formal- 
ism (Spruit 1999, 2002), is still debated: recent works have 
shed some doubts on the efficiency of the Tayler-Spruit dynamo 
mechanism. For instance, Zahn et al. (2007) have shown that al- 
though the Pitts and Tayler instability indeed develop, it does 
not succeed in driving a dynamo. Moreover, this dynamo needs 
a pre-existing magnetic-field seed, and it is not yet clear whether 
magnetic fields were present in the early Universe or not (see 
Giovannini 2004, for a detailed discussion on this subject). 

We presented pre-supernova yields for our stellar models. 
We included them in a chemical evolution model of the Galactic 
halo. Our results predict that the N/0 and the C/O ratios should 
drop at extremely low metallicity (Z < 10"*). For the time being 
though, the observational discrimination between the two curves 
is far beyond our observational possibilities. We also gave the ro- 
tational characteristics of the iron core of our models and showed 
that since they do not lose much angular momentum during their 
evolution, they keep high rotation rates in the core, so their fi- 
nal explosion may be much more energetic than sometimes as- 
sumed. We showed that the most massive of our models may 
contribute to the emission of the gravitational wave background. 

Many interesting questions and issues were not treated in 
this study. Our mass coverage was too coarse to constrain ac- 
curately the limiting masses for undergoing a CNO shell boost, 
for example, or for the different fates at the end of the evolution. 
Even without knowing if magnetic fields were indeed present 
in the early Universe, it would be interesting to study the role 
they could play in primordial stellar evolution. As mentioned 
above, they could drastically change the picture presented here. 
We have also left for a forthcoming paper the question of the 
re-ionisation. Rotation changes the luminosity, the Tetf and the 
lifetimes of the models, so it should change the ionising power 
of the first stars. 
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